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The double integral representing the entropy Stri of spanning trees 
on a large triangular lattice is evaluated using two different methods, 
one algebraic and one graphical. Both methods lead to the same result 
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1 Introduction 



It is well-known that spanning trees on an n-vertex graph G can be enu- 
merated by computing the eigenvalues of the Laplacian matrix Q associated 
with G as 

n— 1 



n$t{g) = - n A 



n 
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where A*, % = 1, 2, • • • , n — 1, are the n — 1 nonzero eigenvalues of Q pQ. 
In physics one often deals with lattices. For large lattices the number of 
spanning trees N$t grows exponentially in n. This permits one to define the 
entropy of spanning trees for a class of lattices G as the limiting value 

S G = lim — In Nst(G). (2) 



The resulting finite number Sg is a quantity of physical interest. 

For regular lattices in d dimensions the eigenvalues A* can be computed 
by standard means |2] leading to expressions of Sg in terms of <i-dimensional 
definite integrals. The resulting integrals, which are independent of boundary 
conditions, have been given for a number of regular lattices [3] , but very few 
of the definite integrals have been evaluated in closed forms. 

To be sure, the entropy for the square lattice has been known and com- 
puted [I| . It is found to be 
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is the Clausen's function and G is the Catalan constant. The algebra re- 
ducing the double integral in Q into the series given in the third line is 
straightforward and can be found in |B] . 
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The only other published closed-form numerical result on the entropy of 
spanning trees is that of the triangular lattice [7j for which one has 
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However, no details leading from the double integral in (jSJ) to the final series 
expression have been published. The expression in the second line is new. 

The purpose of this note is to provide some details of the steps leading 
to the second and third lines in (j5J) which, as we shall see, are nontrivial. We 
evaluate the integral in (jSJ) using two different methods: a direct evaluation 
and a graphical approach. We first evaluate the integral algebraically, and 
then in an alternate approach we use graphical considerations to convert the 
spanning tree problem to a Potts model which is in turn related to an F 
model solved by Baxter 8j. In either case, the same series expression in (|HJ) 
is deduced. 1 



2 Algebraic approach 

In this section we evaluate the integral (0) directly. 
Using the integration formula 
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one of the two integrations in (jSJ) can be carried out, giving 
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Next, we use the identity 
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1 The numerical values of Sm reported in and [H] contain typos in the last two digits. 
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where v = sin (0/2). Then (JBJ becomes 

5trf = In Q+Ji + Ja, (8) 

where 

J m = -jT d# ln(msin#+ vWsin 2 ^), m = l,2 (9) 

after a change of variable 0/2 — ► 0. 

To evaluate Ji, we consider the more general integral 
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The reduction of f'(a) is effected by noting that (V sin 2 9 + a 2 + sin#) 1 = 
(Vsin 2 9 + a 2 — sin 9) /a 2 . Hence, f(a) = (2/7r)Ti 2 (a) and 
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is the inverse tangent integral function described in 
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To evaluate 7 2 , we introduce the following identity 



2 sin^ + V3 + sin 2 # = y/ScosO, 
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Substituting (|T3|) into J 2 in (0) and making use of (JTUJ), we obtain 
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after an obvious change of integration variable in the last step. 
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Next we express the factor in the integrand as 
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and set 
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This yields 
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To express I\ and J 2 in terms of the Clausen's series (J3J), we have from 
Eqs. (2.6), (4.31) and (4.18) of the identities 
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Setting y = \f?> in (JT7|) and 6 = n/3 in (fTHj) . we obtain after making use of 
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Thus, combining (fTTj) and (fTB*j) with (jHJ), we obtain 

5 



Ch 



tri — — ^"-2 
7T 



as given in the second line in (0). 

To express S tTi in the form of the series given in (|SJ), we note that 
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Separating terms in S2 with odd and even denominators, we have 

S2 = Si — — S2 ■ 

It follows that S 2 = (4/5) Si hence S = (6/5) Si, and from J23J and (J22j) we 
obtain 
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This is the result given in (JHJ). 

3 Graphical approach 

In this section we evaluate the integral (0) with the help of a graphical 
mapping. 

The medial graph (or the surrounding lattice) £' of a triangular lattice 
jC is the kagome lattice shown in Fig. 1, where we have shaded its triangular 
faces for convenience. We start from a well-known equivalence between the 
g-state Potts model on £ and an ice-rule vertex model on £' [HI] - [l"2"] - 
Denote the respective partition functions by Zp otts and Zkag-ice- Then the 
equivalence reads 

Zpotts(<?, y/q x) = g n/2 Z kag _ ice (l, 1, x, x, r 1 + xt 2 , t + xt~ 2 ) (23) 

where 

x = (e K - , 
t 3 + t~ 3 = y/q, (24) 

K being the nearest-neighbor interaction of the Potts model, and the ar- 
guments of Zk a g_i C e are weights 2 of the ice-rule model as depicted in Fig. 
2. 



2 In comparing vertex weights with those in it should be noted that for the tri- 
angular lattice which we are considering, roles of shaded and unshaded faces are reversed 
from those indicated in 1121. 
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On the other hand, it is also known that by taking an appropriate q — > 
limit the Potts partition function generates spanning trees [T3]. Write 
the spanning tree generating function in a slightly more general form by 
attaching a weight x to each line segment of the tree. Since the number of 
lines in a spanning tree on an n- vertex graph is fixed at n — 1, this merely 
adds an overall factor x n ~ x . Then, the equivalence reads [7| IT3] 

a^JVsT = limg-(" +1 )/ 2 Z Potts (g, yjq x). (25) 

The equivalences ()23|) and (j23|) hold quite generally for finite lattices, pro- 
vided that appropriate vertex weights are defined for the boundary [11J. 

We next combine (}2*3*j) and (|2*3j) and take the thermodynamic limit n — > 
oo. Assuming the two limits on n and q commute, this leads to 

Stti = -lna; + z kag _i ce (l,l,x,x,t" 1 + xt 2 ,t + xt~ 2 ), t = e in/6 (26) 

where 

ice lim Ti In ^kag— ice • 

We note in passing that the vertex weights in (J26j) now satisfy the free-fermion 
condition [Hj 

1 ■ 1 + x ■ x = (t _1 + xt 2 ) ■ (t + xt~ 2 ) 

so Zkag-ice can be evaluated by standard means such as using Pfaffians [Jj. 
This verifies the integral expression of Stri given in (jHJ. 

Alternately, using a mapping introduced by Lin l-"> (see also [IE1), the 
kagome ice- rule model can be mapped into an ice-rule model on the triangular 
lattice. The trick is to 'shrink' each shaded face of the kagome lattice into a 
point as shown in Fig. 1. This results in a 20- vertex ice- rule model on the 
triangular lattice which has 3 arrows in and 3 arrows out at every vertex. 

There are 4 different vertex types in the 20-vertex model: (i) six vertices 
in which the 3 incoming arrows are adjacent, (ii) two vertices in which the 
incoming and outgoing arrows alternate around the vertex, (iii) six vertices 
in which 2 incoming arrows are adjacent while the third in- arrow points 
in a direction opposite to one of the 2 adjacent incoming arrows, and (iv) 
six vertices in which 2 incoming arrows are adjacent while the third in-arrow 



8 



points opposite to the other adjacent incoming arrow. These 4 configurations 
are shown in the left-side column in Fig. 3 to which we assign respective 
energies e», i = 1,2,3,4. The other 16 ice-rule configurations are obtained 
by rotations. 

The shrinking processes from which these ice-rule configurations are de- 
duced are also shown in Fig. 3. Write the 4 weights as a, b, c, c', respectively. 
We read off from Fig. 3 to obtain 
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where Ztri-ice is defined similar to Zkag-ice- 

Baxter jH] has considered the evaluation of Ztri-ice in the subspace of c = c'. 
Now the right-hand side of (J2l)|) is independent of the value of x so we have 
the freedom to choose a value of x for which z ice _ tr i can be evaluated. It is 
readily verified that by taking x = l/y/3 and t = e J7f//6 we have 

a=l/\/3, b = 3a, c = c = 2a 

which is a case solved by Baxter. In this case Baxter obtained 
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where P denotes the principal value. Combining (J26j) and (J28|) with (J29j) and 
a = l/\/3, we obtain 
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This is the result given in (j5J) as first reported in jjj. 
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Figure captions 

Fig. 1. The shrinking of a kagome lattice into a triangular lattice. Each 
shaded triangular face is shrunk into a point. 

Fig. 2. The six ice-rule configurations of the kagome lattice and the associ- 
ated vertex weights. 

Fig. 3. The four types of ice-rule configurations of the triangular lattice and 
the associated shrinking processes. 
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